Effective forces in colloidal mixtures: from depletion attraction to accumulation 

repulsion 
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Computer simulations and theory are used to systematically investigate how the effective force 
between two big colloidal spheres in a sea of small spheres depends on the basic (big-small and 
small-small) interactions. The latter are modeled as hard-core pair potentials with a Yukawa tail 
which can be both repulsive or attractive. For a repulsive small-small interaction, the effective force 
follows the trends as predicted by a mapping onto an effective non-additive hard-core mixture: both 
a depletion attraction and an accumulation repulsion caused by small spheres adsorbing onto the 
big ones can be obtained depending on the sign of the big-small interaction. For repulsive big- 
small interactions, the effect of adding a small-small attraction also follows the trends predicted 
by the mapping. But a more subtle "repulsion through attraction" effect arises when both big- 
small and small-small attractions occur: upon increasing the strength of the small-small interaction, 
the effective potential becomes more repulsive. We have further tested several theoretical methods 
against our computer simulations: The superposition approximation works best for an added big- 
small repulsion, and breaks down for a strong big-small attraction, while density functional theory 
is very accurate for any big-small interaction when the small particles are pure hard-spheres. The 
theoretical methods perform most poorly for small-small attractions. 
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I. INTRODUCTION 



A fundamental description of colloidal interactions 
based on statistical mechanics is needed to understand 
and predict the stability and the phase behavior of col- 
loidal suspensions Most inter-colloidal forces are 
effective in the sense that some microscopic degrees of 
freedom are averaged out. This concept of mean or ef- 
fective interactions [f|-[| is crucial to bridge the different 
length scales involved in colloidal systems and has been 
exploited in many different circumstances. Examples of 
the microscopic degrees of freedom include solvent parti- 
clesjl , smaller colloidal particles , added polymer coils 
|^| |llj or monomers of (grafted) polymer chains 12-]1| 
as well as counter-ions [ p^|Jl7| ] and salt ions [Q in the 
case of charged suspensions. The resulting effective forces 
turn out to exhibit a wide range of features: They can 
be attractive, repulsive, or oscillatory, and are an impor- 
tant key to understanding colloidal stability as well as 
flocculation and coagulation [EJ. 

In the present paper we consider a binary colloidal 
mixture of big and small colloidal particles and inves- 
tigate the distance-resolved effective force and potential 
between two big colloidal spheres surrounded by many 
small colloidal spheres, as depicted in Fig. [l]. The small 
particles are averaged out so that one is left with an ef- 
fective pair potential between the big ones, which is su- 
perimposed onto the direct big-big interaction. The form 
of these effective interactions depends sensitively on the 
basic big-small and small-small interactions. 



Most previous studies have focused on a hard-sphere 
(HS)-like interaction between the colloidal particles, with 
a special emphasis on two cases: (1) the Asakura-Oosawa 
model, which was originally designed to describe mix- 
tures of sterically-stabilized colloids and nonadsorbing 
polymer coils ||. In this strongly non-additive |ls| ] 
model, the big-small interaction is HS like, while the 
small-small interaction is ideal gas like. The effective pair 
potential between two colloids can be exactly calculated 
and has an attractive well proportional to the density of 
the small particles, and a range equal to their diameter. 
(2) Much work has also focused on the complementary 
model of additive HS mixtures, inspired in part by a pre- 
diction that they might phase-separate pQ]. The effec- 
tive depletion forces were obtained by computer simula- 
tion pUJ^I , approximate theories | |23||24f| , density func- 
tional theory (DFT) Hill, and experiments ff|7|||§. 
Simplified potentials were used to investigate the phase- 
behavior of binary HS mixtures [ p9| , where it was found 
that fluid- fluid phase-separation, when it occurs pfj| , ^l| , 
was always metastable w.r.t. fluid-solid phase-separation. 

In this work we generalize these studies to arbitrary 
big-small and small-small interactions. Our study is mo- 
tivated by the fact that real colloidal suspensions typ- 
ically have soft interactions which are beyond the HS 
model. These soft interactions can easily be tuned, for ex- 
ample by changing the solvent composition, the amount 
of added salt, the surface charge, etc., which provides a 
route to tailoring the effective interactions between the 
big particles. Hence it is essential to understand on a full 
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statistical level how the basic big-small and small-small 
interactions affect the effective big-big interactions in a 
binary colloidal mixture. 

There has been some previous work on how non-HS 
like interactions affect the effective big-big interactions, 
see e.g. refs. |3,^2|-^3|, which showed that the differences 
with the pure HS case could be substantial. However, 
with a few exceptions, all these studies were performed 
with approximate methods such as integral equations. 
Our goal here is to provide benchmark exact computer 
simulation results for a set of 9 different parameters cor- 
responding to combinations of attractive and repulsive 
small-small and big-small interactions. This allows us 
to systematically study the trends in the effective forces 
as induced by the basic interactions. It also allows us 
to test several theoretical techniques, namely DFT, the 
superposition approximation 44 1, and a mapping to a 
non-additive HS model ]3C| ]. 

Generalizing from an additive to a non-additive HS 
mixture already leads to a much richer class of effective 
potentials ^Q.^O 31 1 , which can be exactly calculated 
from an accurate theory of additive HS mixtures [ p6[ . 
Similarly, we find here that the effective big-big inter- 
actions are profoundly influenced by the basic big-small 
and small-small interactions. The trends can be sum- 
marized as follows: Adding a big-small repulsion results 
in enhanced attraction through the standard depletion 
mechanism J|||. Adding a big-small attraction leads to 
an accumulation of the small particles near each big one, 
which in turn results in a more repulsive big-big inter- 
action (accumulation repulsion). For a given big-small 
interaction, adding a small-small repulsion also results 
in an enhanced small particle density near the big parti- 
cles, and therefore in a more repulsive big-big interaction. 
Adding a small-small attraction when the big-small in- 
teraction is repulsive results in a more attractive effective 
interaction. All these effects can be qualitatively under- 
stood from a mapping to a non-additive HS Model. But 
when a small-small and a big-small attraction are com- 
bined, a repulsion through attraction effect occurs, which 
is not captured by DFT or by our mapping scheme. 

The paper is organized as follows: In Chapter II, we 
discuss how to formally map a two-component mixture 
onto an effective one-component one, and define our tar- 
get quantities. In Chapter III, we describe our model for 
the basic interactions. We discuss our simulation results 
in Chapter IV, and the results of several different theo- 
ries in Chapter V. Our conclusions are stated in Chapter 
VI. 



II. MAPPING A BINARY MIXTURE ONTO AN 
EFFECTIVE ONE-COMPONENT SYSTEM 

In this section we briefly describe how to map a 
two-component binary mixture onto an effective one- 
component system. This statistical mechanical proce- 
dure was originally developed by McMillan and Mayer 



|47[ to trace out the solvent in a suspension, but it can 
be used just as well to trace out a smaller component in 
a binary colloid mixture. 

We begin with the Hamiltonian of our two-component 
mixture 



H = K + H i 



bb 



bs ■ 



(1) 



Here K is the total kinetic energy of the mixture and 
three potential energy contributions are: 



JV ( , 



i<j 

if ss = X> ss (rf-rj), 



i=i i=i 



(2) 



where the $y denote the pairwise interaction potentials, 
and v b i and r s 4 denote the coordinates of the centers of 
the big and small particles respectively. 

The formal procedure to reduce the two-component 
system to an effective one-component one, described in 
more detail elsewhere, see e.g. |4p|. p^| , p9| , ^l| , consists of 
tracing out the small particles for a fixed configuration 
{r^} of the big particles. This is most conveniently done 
for a semi-grand ensemble where the small particles in 
the volume V are kept at a fixed chemical potential fx s . 
The resulting effective one-component Hamiltonian that 
governs the behavior of the big particles takes the form: 



Htl S = H bb 



(3) 



where fl = Q(Nb, z 3) V; {r^}) is the grand potential of 
a fluid of small particles at fugacity z 8 = Xj 3 exp(/3/x s ), 
(/3 = l/ksT) subjected to the external potential of the 
fixed configuration {r^} of the big particles. This grand- 
potential can be further expanded as a sum of (n-body) 
terms |tJ9|: 



n=0 



(4) 



The first term, Oo, is the so-called volume term, which 
only depends on properties of the small-particles (i.e. 
$ ss ) and is therefore independent of the configuration 
{r^}. The one-body term fii can be related to Nb times 
the free-energy gained by inserting a single big particle 
into a sea of small particles. It depends on both Q ss and 
& bs , but is independent of $ bb - Of most interest to us 
here is the two-body term f^, which can be written as: 



n 2 (N b ,z s ,{rt}) = ^V bb 

Kj 



eff 



(5) 
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where the two-body effective pair-potential V£f^(r) is 
related to the free-energy difference between two big 
spheres a distance r apart, immersed in a sea of small 
spheres at fugacity z s , and the same system with the 
two spheres at r — oo. It is important to notice that 
V h \ (r) itself only depends on $f, s and <& ss , and not 
directly on This observation was a key to under- 

standing why the effective potentials for non-additive HS 
mixtures can be calculated by a theory for additive HS 
mixtures @|l|,|6|. Nevertheless, through Eq. (|), the 
effect of $bf,(r) is felt in that it sets the range where 
V bh (r) is relevant. 

In a similar way higher order interactions can be de- 
rived ]lq ], but for this scheme to remain tractable, one 
usually truncates at the pair level. As long as the ranges 
of $ ss and $b s are small compared to the range of $(,6) 
higher order terms are not expected to be very important 

Thus far we have developed this formal tracing out in 
the semi-grand ensemble, which helps emphasize that in 
a mixture with many large particles, the correct effective 
potential V b e J^ (r) to be used in Eq. (||) is the one fixed 
by the chemical potential or fugacity instead of the over- 
all density of the smaller particles, as first emphasized 
by Lekkerkerker et al. jn|. It is important to keep this 
in mind when studying phase-behavior. However, for a 
practical calculation of the pair term, with only two big 
particles in an infinite sea of small ones, one can just as 
easily work in the canonical ensemble, keeping in mind 
that when using the obtained V bb (r) for a system with 
many big particles, the input density p s is really that of 
a reservoir of small particles kept at the same chemical 
otcntial as the full mixture of big and small particles 
!§§■ 

For the case of two big-spheres a distance r apart, in 
a bath of small particles at density p s , the average force 
is related to the effective potential by: 



bb 



(r) 



d_ 
Or 



V b l ff (r)- 



(6) 



This average force provides an intuitive way of under- 
standing the effective interactions which parallels the 
more formal derivation above. First we define a one-body 
density as follows: 



r') =«£«5(r'-r?)», 



(7) 



where <<>> denotes a canonical thermodynamic aver- 
age over the small particles. If one big particle is fixed 
at the origin, and one is fixed at r, then we can define 
a one-particle density of the small particles pi 1 ^ (r'; 0, r) 
for that fixed configuration. The effective force induced 
between two big spheres by the smaller particles can then 
be written as |44|]21[ |: 



bb 



o,r0 



d_ 
dr' 



<kb s {r')dv' ; 



(8) 



which has the obvious physical interpretation that the 
total effective or mean force is simply the average of the 
sum over all the big-small interactions. 

In the next sections, we will concentrate on the effect 
of changing $ ss (r) and $b s (r) , at fixed reservoir density 
p s , on the effective interactions Fu/^(r) and (r). 



III. MODELS 

A. Interactions 

We model the basic small-small and big-small interac- 
tions defined in Eq. (||) as hard-core Yukawa pair poten- 
tials: 



*«(r) 



r < a. 



<f>ij (r) otherwise, 



(9) 



where in each case r denotes the distance between the 
centers of the relevant particles. Throughout, we take 
4>bb{r) = 0, while the big-small interaction added to the 
HS repulsion is given by: 



4>bs(r) = bs(Tbs exp [-Kb s (r - a bs )} 
r 



(10) 



and the small-small interaction added to the HS repulsion 
is given by: 



cj> ss {r) 



■ exp [-n ss (r - a ss )} 



(11) 



where a bs = \{a bb + a ss ). 

The Yukawa tail can be either repulsive or attractive. 
This gives us the possibility to change the interaction 
over a broad range from pure HS (which serve as a ref- 
erence case) to soft repulsions and attractive tails. We 
note that a hard-core Yukawa interaction has been fre- 
quently used in theoretical studies on liquids, see e.g. pTjj ] 
and refs therein. Both an attractive and repulsive tail 
is realized in various colloidal solutions: a short-ranged 
attractive Yukawa tail has been shown to satisfactorily 
model the "stickiness" in the interaction between globu- 
lar protein solutions |50|. A repulsive Yukawa tail, on the 
other hand, describes charged suspensions |§,|lj where 
the range is controlled by the added salt concentration 
and the amplitude by the colloidal charge. 



B. Parameter combinations 

Our aim is to systematically investigate the effect of 
the basic big-small and small-small interactions on the 
effective big-big interactions. We vary the big-small in- 
teraction by changing 0e bs , and vary the small-small in- 
teraction by changing (3e ss . The other parameters are 
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fixed as follows: a ss = 0.2t7b6, i^bs — 6/cbb — 1-2/<t ss , 
and k ss — 15/crbb = 3/cr ss . The range of the big-small 
interaction is of the order of the small particle size, while 
the range of the small-small interactions is significantly 
less than the small particle size. For all simulations the 
packing fraction of the small particles in the bulk is set 
to T) s = irp s a 3 s /6 = 0.1. 

In total, we studied 9 different parameter combinations 
corresponding to taking f3eb s and f3e ss to be positive, neg- 
ative, or zero. The detailed values are given in table [j]. 



IV. COMPUTER SIMULATIONS OF EFFECTIVE 
FORCES AND INTERACTIONS 

A. Simulation method 

Our simulation set-up contains two big spherical par- 
ticles in a large cubic simulation box of length L with 
periodic boundary conditions in all three directions. The 
big particles with their center-to-center separation r are 
fixed along the body diagonal of the cubic box and 
N s = 10000 small mobile particles of diameter <j ss are 
added to the box. The box length L — 7A9a ss is suffi- 
ciently large to exclude any spurious periodic image ef- 
fects. We also studied a single big particle in the sim- 
ulation box to access the one-body small-particle den- 
sity p s { r ) around an isolated big sphere. For two big 
spheres p s — N S /(L 3 — ^Trcr 3 b ), while for one big sphere 
p s = N S /(L 3 — ^TT(J^ b ). We checked the asymptotic den- 
sity profiles to confirm that the bulk densities were always 
the same. 

We use a molecular dynamics (MD) code combining 
the velocity Verlet algorithm S| with discrete collisions 
and reflections induced by the hard core of the interac- 
tion potentials in order to calculate the trajectories of the 
small particles. The system was carefully equilibrated 
and then statistical averages were computed such as the 
mean force acting on the big particles (see Eq. (Q)) or the 
inhomogeneous density field p s (r) of the small particles 
around a single big particle. Details of our simulation 
procedure are described in references |lg|j53|] . 

B. One body density profiles 

Before we discuss the depletion potentials it is useful 
to first examine the density profiles p s (r) of the small 
particles around a single sphere. 

1. No extra big-small interaction: runs 1-3 

In Fig. ^ these are shown for runs 1 — 3, where ebs = 0. 
Compared to the density profile of the pure HS, there is a 
significantly increased accumulation near the big spheres 
when a small-small repulsive interaction is added. To 



first order this can be understood by mapping the <p sa (r) 
onto an effective HS diameter a e JJ Q: Adding a small- 
small repulsive interaction increases the effective HS size, 
and therefore also the effective HS packing fraction, re- 
sulting in a more pronounced accumulation of density 
near the surface. Adding a small-small attraction has 
the opposite effect. This is because the bulk is now more 
favorable for the small particles, and they are attracted 
to it. 



2. Added big-small repulsion: runs 4~6 

In Fig. the density profiles are shown for runs 4 — 6, 
where (3tbs = 0.82. As expected, the repulsive big-small 
interaction leads to a reduction of the density at contact. 
Again, adding a small-small repulsion increases the den- 
sity at contact, and adding a small-small attraction leads 
to a further depletion of the density at contact. 



3. Added big-small attraction: runs 7-9 

In Fig. |], the density profiles are compared for runs 
7 — 9, where fle^s = —0.82. As expected, the attrac- 
tive big-small interaction results in an enhanced density 
at contact. Similarly to the previous two cases where 
€b s = and ebs > 0, respectively, adding a small-small re- 
pulsion enhances the contact density w.r.t. the case of no 
small-small repulsion. But in contrast to the two previ- 
ous cases, where adding a small-small attraction resulted 
in a depleted density profile w.r.t. the pure HS case, here 
adding e ss < results in an enhanced total accumulation 
of the small particles near the big one. Although the con- 
tact value is slightly smaller than the case for pure HS, 
there is a marked accumulation further out, correspond- 
ing to a second layer of particles. The relative adsorption 
of the small particles around the big one is thus consid- 
erably larger than for the case of no small-small attrac- 
tions. This can be understood by the following simple 
argument: the big-small attraction leads to an accumu- 
lation of the small particles near the surface of the big 
sphere. When this accumulation is large enough, it be- 
comes favorable for the small particles to leave the bulk, 
and approach the surface of the big particle, where their 
local density is larger. Thus the two attractions amplify 
each other in a non-linear fashion. 

In conclusion then, when comparing Figs. [^, || and ^, 
it is clear that for the parameters we have chosen, the 
big-small interaction has the largest relative effect on the 
density profiles. The effect of adding a small-small repul- 
sion can be qualitatively understood by the larger effec- 
tive sphere size and concomitant larger packing fraction. 
The effect of adding a small-small attraction can be qual- 
itatively understood by the fact that the bulk is usually 
preferred over the surface of the particle, except when 
the big-small interaction is strong enough to provoke a 
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non-linear enhancement of the density of small particles 
near a big one. 

Having investigated the effect of the interactions on 
the one-body density profiles, we now turn to the related 
two-body depletion potentials. 

C. Effective pair forces and potentials 

In Fig. H we compare the effective depletion force be- 
tween two big spheres for all nine parameter combina- 
tions detailed in table |[ In Fig. ^ we compare the re- 
lated depletion potentials. First we note that changing 
the big-small interaction is the dominant effect: the de- 
pletion forces and potentials split naturally into three 
groups: no added big-small interaction (solid lines, runs 
1 — 3), big-small attraction (dotted lines, runs 4 — 6), and 
big-small repulsion (dot-dashed lines, runs 7 — 9). We 
will treat each case in turn. 



fact that the big-small attraction results in an enhanced 
density of small particles near a single large sphere. 
When two such large spheres approach one another, the 
layers of small particles around each one begin to inter- 
act, leading to an enhanced repulsion between them. As 
demonstrated in Fig. ^ for a given big-small attraction, 
adding small-small repulsion or attraction both result in 
a further enhanced density p s (r) of the small particles 
around a single large one, which in turn explains why 
the effective pair potentials are more repulsive for both 
e ss > and e ss < 0. Relatively speaking, adding small- 
small attractions has the largest effect on the effective 
potentials, which is due to the non- linear coupling be- 
tween the big-small and the small-small attractions. 

V. THEORETICAL DESCRIPTIONS OF THE 
EFFECTIVE FORCES AND POTENTIALS 

A. Superposition approximation 



1. No extra big-small interaction: runs 1-3 

As was already seen for the one-body profiles, adding a 
small-small repulsion (run 2) results in a larger effective 
small-sphere size and packing fraction, which is reflected 
in more pronounced oscillations compared to the pure HS 
case (run I). These are evident both in the effective force 
and in the effective pair potential. Adding a small-small 
attraction (run 3) results in a reduced density near a big 
particle, as seen in Fig. ||. This would imply that each big 
sphere excludes slightly more free volume than if there 
were no small-small attractions, which, in turn, implies 
a slightly more attractive potential, as is observed. 



2. Added big-small repulsion: runs 4~6 

The dominant effect of adding a big-small repulsion is 
to make the effective forces and potentials much more 
attractive. This can be understood with the classical 
picture of depletion H,||: Adding a big-small repulsion 
results in a larger depletion layer around each large par- 
ticle. When two large spheres approach, the amount of 
doubly excluded volume is therefore larger, resulting in 
a more attractive effective (depletion) potential or force. 
Again, adding a small-small repulsion (run 5) results in 
enhanced layering as compared to the pure HS small par- 
ticles. Adding a small-small attraction (run 6) has only 
a weak effect similar to what was seen for run 3. 



3. Added big-small attraction: runs 7-9 

The dominant effect of adding a big-small attraction 
is to make the effective forces and potentials much more 
repulsive. This can be qualitatively understood from the 



To calculate the effective forces via Eq. (g), one needs 
a prescription for calculating the one-body density of the 
small particles. In section IV, we essentially did this by 
computer simulations. In this section we approximate the 
full one-body density by a superposition of the one-body 
density p s (r) around an isolated single sphere E3): 



pW(r';0,r) = p s (r')p s (\r-v'\)/p s , 



(12) 



an approach similar in spirit to the Kirkwood superpo- 
sition approximation (54). Since all input information 
comes from the (radially symmetric) problem of a single 
sphere, this superposition approximation greatly simpli- 
fies the calculation of the two-body depletion forces. The 
input p s (r) could come from density functional theory or 
integral equation theory, as was done previously by other 
authors |^l],[34j], but here we will use the Psir) generated 
by our simulations and depicted in Figs. prfl- The results 
are shown in Figs. [^15|, where we compare in detail the 
radial force F^^(r) — F^jf (r) • r/r and the effective po- 
tential VSf (r) calculated with the superposition approx- 
imation, to the results obtained by direct simulations. 

Figure [7] (run 1) shows the HS reference case. Here 
the superposition approximation works quite well. The 
packing fraction of the small spheres rj s =0.1 is rather 
low, so we expect that the small-small correlation effects 
are not very strong. The total one-body density is there- 
fore well approxi mated by Eq. (|I2|). As demonstrated 



by other authors [|21 34 for the pure HS case, this su- 



perposition approximation begins to break down as rj s 
increases and two-body correlation effects become more 
important. For example, they found that for large rj s , the 
superposition approximation leads to an overestimate of 
the strength of the attractive force at contact. This is 
exactly what is seen in Fig. || (run 2), where the effect 
of an added repulsive <f> ss (r) can be understood in terms 
of a larger effective packing fraction 77,5. Figure || (run 3) 
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shows that for e ss < the superposition approximation 
shows a similar error to what was seen for run 2, i.e. the 
forces and potentials are too attractive. 

In Figs. io|-|l2] (runs 4-6), where the repulsive eb s (r) 
induces much more attractive effective interactions than 
for the pure HS case, the superposition approximation is 
seen to work quite well for e ss = and e ss > 0. This 
is most likely because the effective (depletion) force or 
potential is dominated by 4>bs{f), which also induces a 
lower p s (r) (see Fig. ||), so that small-small correlations 
only play a relatively minor role and Eq. ( p"2| ) is rather 
accurate. Even for e ss < 0, the results are quite good, 
although the superposition approximation tends to pre- 
dict forces and potentials that are too attractive, just as 
was found for 6bs{r) = 0. 

In Figs. [13|-[l5| (runs 7-9), where the attractive 0;, s (r) 
induces an increased local accumulation of the small 
spheres near each big sphere, the superposition approxi- 
mation performs less well. This is perhaps not surprising, 
as the enhanced density of small particles near each big 
sphere results in a more important role for small-small 
correlations, which are not well treated by the superpo- 
sition approximation. 

In summary then, for HS or other repulsive small-small 
interactions, the superposition approximation works best 
for low (effective) packing fractions rj s , where correlations 
between the small particles do not strongly alter the full 
one-body density p^P{r ! ; 0, r) from the superposition of 
the densities around an isolated big sphere. The case 
of a strongly repulsive 4>b s (r) , which lowers the effective 
packing fraction near the big spheres, is particularly well 
described by the superposition approximation. In con- 
trast, the case of an attractive 4>bs{r), which results in 
an increased accumulation of small particles around each 
big sphere, is not as well described. 



B. Density Functional Theory 

In a recent development, density functional theory 
(DFT) has been used to derive effective potentials for 
additive |25@ and non-additive |o|j3|,||l HS mixtures. 
For the additive case, quantitative agreement with com- 
puter simulations was achieved |2^ , |2q| . Since the calcu- 
lations for the non-additive case were shown to be equiv- 
alent to the additive ones fl30| , [3l"l| , a similar quantitative 
accuracy could be expected there. In brief, the method 
is based on the following exact relationship between the 
effective pot ential and the one-body direct correlation 
function [E5l,E6|: 



c«(r)=-/3 



(3V b % ff {v)= lim (c^(oc) 

Mb— > — CO 



cl 1} (r)) 



(13) 



where (r) is defined for the case where a big sphere is 
fixed at the origin, and exerts a field on the small spheres 
and on a big test particle inserted at r Q. DFT can 
therefore provide a route to the effective potentials since 



SFex [Pb, Ps 

Spb{r) 



(14) 



where ^F ex [PbiPs] is the excess (over ideal) intrinsic 
Helmholtz free energy functional of the mixture [ pBI . 
Given some mixture functional, one can obtain the ef- 
fective potentials and forces from a radially symmetric 
calculation of c£ ( r ) m the limit of vanishing density of 
the big spheres. This is much simpler than say trying to 
use DFT to directly calculate the cylindrically symmet- 
ric one-body density p^p (r'; 0, r), for use in Eq. (||) |^6| . 
All one now needs is some prescription for the mixture 
functional. Here we use Rosenfeld's very successful fun- 
damental measure theory |5(|, in its form valid for HS 
mixtures. As found previously |2q] , the DFT and sim- 
ulation results for effective interactions of the pure HS 
case depicted in Fig. ^ (run 1) agree very well. 

This DFT method can easily be extended to an arbi- 
trary cj)b s (r), since this simply corresponds to an addi- 
tional external field in Eq. ([13]). Previous DFT calcula- 
tions f30|| found good agreement with earlier simulations 
with an attractive <fib s (f) ]35| . Here we also find very good 
agreement with the runs 4, and 7, (Figs. |l^ and [l3|) which 
correspond to a finite 4>b s (j") but no additional </) ss (r). 

Since there is at present no successful two-component 
DFT for mixtures where the small-particles are not HS 
like, some approximations must be made. We chose to 
map the ss (r) onto effective HS diameters using 
the Barker- Henderson approach [b7| : 



eff 



(1 - exp [-/%,- (r)])dr, 



(15) 



where cr^ is the bare HS diameter. We then calculated 
effective pair potentials and effective pair forces with our 
full DFT approach, including an explicit 4>bs{r) when 
needed, but with the small-spheres mapped onto the ef- 
fective diameters a^J = 1.447(T SS and — 0.658er ss 
for the repulsive and attractive ss (r) respectively. Simi- 
larly the effective packing fraction becomes rfJJ = 0.303 
for the repulsive and rfJJ = 0.0290 for the attractive 
interaction. 

Firstly for runs 2, 5, and 8 (Figs. |, [ll], and[l4[ respec- 
tively), which all correspond to a repulsive </> ss (r), we find 
good agreement for no additional <pb s {r) (see Fig. ||) , but 
less good agreement for a repulsive or attractive <fib s (,r) 
(see Fig. [ll] and Fig. [l4| respectively) . 

Because of the strength of the small-small attraction, 
the mapping results in a very low effective packing frac- 
tion of the small particles. Overall, the DFT underes- 
timates the effective forces for tb s — and e& s > 0, 
as can be seen in Figs. ^| and [l^. It performs rather 
poorly for run 9, where the attractive 4>bs(r) results in 
a non-linear enhancement of the small-particle density 
profile, an effect not taken into account with our HS 
mapping. This suggests that a different two-component 
DFT, which explicitly takes into account the small-small 



G 



attraction ifissir) needs to be developed before we can 
use this route to derive accurate effective pair forces and 
potentials. 



C. Mapping to non-additive HS system 

In a previous paper |]30| two of us proposed that the 
effects of big-small and small-small interactions on the 
V b e b (r) could be understood by mapping onto those of 
non-additive HS systems. These have the advantage that 
they can be determined by an exact mapping onto the de- 
pletion potentials of additive HS mixtures ]3(^,[3l| , which, 
in turn, are well understood and for which a good param- 
eterization exists ^6|. Even a small non-additivity was 
shown to have a large effect on the depletion potentials. 

By mapping the big-small and small-small interactions 
onto effective HS diameters through Eq. ([l5]) one can de- 
fine the non-additivity in terms of the parameter A: 

af s f ^^(a bb + afJ)(l + A). (16) 

Four different ways of adding interactions to introduce 
non-additivity were studied [p0| : 

(i) repulsive (3cf) ss (r): a e b f J = a bs ; a e JJ > a ss ; A < 

(ii) attractive f3(j) ss (r): a e JJ = a bs ; a e JJ < <r ss ; A > 

(iii) repulsive f3(j) bs (r): a^ f > a bs ; a e JJ = a ss ; A > 

(iv) attractive /3<t> bs (r): u e JJ < a bs ; a e JJ = a ss ; A < 0. 
For details of how each of the different cases of non- 
additivity affect the depletion potentials, we refer to refs. 

©ma- 

An example where the mapping to non-additivity 
works well is given in Fig. |^ (run 2), where we mapped 
the small-small repulsion onto an effective HS diameter 
(case (i)), so that the DFT result is really that of a non- 
additive HS mixture with A = —0.069. In rcf. |Q we also 
found semi-quantitative agreement with the mapping for 
weak and short-ranged (3(j) bs {r). In the present case, 
where (3(j> ba (r) is stronger and longer ranged, the agree- 
ment is no longer quantitative. Nevertheless, for runs 1-8 
the mapping scheme provides a qualitative explanation 
of the trends. I.e. for a fixed 4> bs (r), adding a repulsive 
(case (i)) or attractive (case (ii)) (f> ss (r) causes V bb {r) to 
become more repulsive or attractive, respectively. Simi- 
larly, for a <f) ss {r), adding a repulsive (case (iii)) or at- 
tractive (case (iv)) (j> bs (r) results in a more attractive or 
repulsive V bb {r), respectively. The only case where this 
scheme breaks down is run 9, where adding </> ss (r) < 
for an attractive 4> bs (r) should fall under case (ii). But 
instead of inducing more attraction, the effective big-big 
potential becomes more repulsive. Of course this is not 
surprising, since the DFT results already showed that a 
mapping scheme misses the non-linear enhancement of 
the small particle density profiles. 

In conclusion then, the mapping to non-additivity 
works best for a repulsive small-small interaction. For 



repulsive and attractive big-small interactions, the map- 
ping is only qualitative. The much better semi- 
quantitative agreement found in ref. |30| can be traced 
to the much weaker effect of the 4> bs (r) used there. Just 
as was seen for the direct DFT methods, it is the case of 
small-small attraction combined with big-small attrac- 
tion which seems most difficult to capture within our 
mapping scheme. 



VI. CONCLUSIONS 

In conclusion, we have shown how the basic interac- 
tions in a colloidal mixture control the resulting effective 
interactions between the big particles. This knowledge 
may be exploited to stabilize colloidal particles against 
coagulation and to tailor the colloidal phase diagrams. 

Adding a repulsive <j) bs {r) results in a strongly en- 
hanced attraction through the standard depletion mech- 
anism. We also found at least two ways to obtain sig- 
nificantly more repulsive effective interactions caused by 
accumulation of the small particles (accumulation repul- 
sion): (1) Adding repulsions between the small particles 
results in an enhanced accumulation near the surface of 
the large particles; when two large particles approach 
each other, this results in an effective repulsion between 
them. (2) Adding an attraction between the large and 
the small particles also results in an enhanced accumula- 
tion near the surface of the large particles and therefore 
in repulsive effective interactions. Furthermore, we found 
that for an attractive big-small interaction, adding small- 
small attractions resulted in even more effective repul- 
sion. This "repulsion through attraction" effect is caused 
by a coupling between 4> ss {r) and 4> b s(r). 

These predictions could in principle be verified in ex- 
periments that directly measure the effective forces of 
colloidal suspensions, such as optical tweezer s Rj5q| o r 
Total Internal Reflection Microscopy (TIRM) J2j,|§|59j. 
They could also be verified indirectly through measure- 
ments of phase-behavior and coagulation. Measurements 
of the second osmotic virial coefficient might also be very 
sensitive probes of the effective interactions pjj ]. Sys- 
tems where these interactions could be tuned include 
for example ternary suspensions |60[| where the small- 
small attraction is generated by depletion attraction of an 
even smaller third colloidal component or highly salted 
charged suspensions where the van-der-Waals attraction 
dominates. 

In some interesting recent experiments fill , a colloidal 
suspension of neutral big particles was stabilized by the 
addition of highly charged small nanoparticles. The pro- 
posed mechanism was termed "nanoparticle halos" , and 
is very similar to the mechanism we observe, for exam- 
ple, in our run 2, where a repulsive 4> S s{r) of the screened 
Coulomb (Yukawa) form was added to the smaller par- 
ticles, resulting in an increased accumulation of small 
particles near each big one, and an effective repulsion 
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between the big particles. The effects in the experiments 
may also be enhanced by a small attractive 4>bs(f) (like 
our run 8). We are currently actively pursuing a more 
detailed comparison with these experiments. 

Another possible application of this work is to super- 
critical solvents |6^] , which have important applications 
in industrial processes. The question of how the effective 
force on the big particles depends on the interactions with 
a low density solvent is encountered there as well [^fll] , 
sometimes for similar size ratios. 

More generally to make useful predictions for practi- 
cal applications, one needs to more thoroughly explore 
the rather large parameter space, which includes the , 
the Kij, rj s and the ratio a ss /<Jbb- This would be very 
tedious with simulations - trustworthy theoretical tech- 
niques would be more practical. We attempted a number 
of theoretical descriptions of the effective interactions. 
The superposition approximation, which has the advan- 
tage of only needing one-body input, works best for a 
strong big-small repulsion, but becomes much less ac- 
curate when the (local) small particle packing is high. 
When the smaller component is purely HS like, then our 
DFT approach is very accurate, just as was found for ad- 
ditive and non-additive HS mixtures (2^j3(],|3l],|^| . When 
a repulsive interaction is added between the small par- 
ticles, DFT does not perform quite as well, and when 
an attraction is added between the small particles, the 
differences are even more important. By mapping the 
small-small and big-small interactions onto effective HS 
diameters, we can map onto an effective non-additive HS 
model. This explains the qualitative trends for most of 
our parameter combinations. It is quantitative if only 
a small-small repulsion is added, but breaks down when 
both a small-small and a big-small interaction couple to- 
gether to induce an enhanced small particle density. 

All three tested theoretical methods perform less well 
for added small-small attractions. Constructing a reli- 
able theory to treat this very interesting case is therefore 
a challenging problem. One possible way to extend the 
DFT calculations to an attractive small-small interaction 
would be to add a mean-field attractive term (55). For 
example, mean-field like functionals based on thermody- 
namic perturbation theory around a HS reference sys- 
tem have been successfully applied to the density profiles 
and phase-behavior in systems with attractive potentials 

New physics is expected when wetting or drying phe- 
nomena control the density of the small particles between 
the big ones. This is relevant if the bulk fluid of the small 
particles is close to liquid-gas phase coexistence. A wet- 
ting transition is expected to have a profound impact on 
the effective interaction as well, one important effect is 
a liquid "bridge" of small particles between the big one 
which has been recently studied in more detail |65f| . 
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TABLE I. Parameter combinations for simulation runs. 
The interactions are given by Eqs. (|lo|) and ([ll]), and only 
the et s and e 3S are changed; the other parameters are kept 
fixed. 



run 


(3e bs 


f3e ss 


1 








2 





2.99 


3 





-0.996 


4 


0.82 





5 


0.82 


2.99 


6 


0.82 


-0.996 


7 


-0.82 





8 


-0.82 


2.99 


9 


-0.82 


-0.996 







0% 




Two big spheres of radius abb experience an ef- 
fective force F^f(r) induced by the sea of small spheres or 
radius a ss ■ 



run 4 e b 


s >o; e B 




a o run 5 e b 


s >o; e s 




□ □ run 6 e b 






run 1 








FIG. 3. Normalized density profiles p a (r)/p s from com- 
puter simulations for runs 4 — 6; ej, s > 0, e sa is varied. We 
also show the result for run 1, the pure HS case. 



4 - 



3 - 



run 1 e bs =0 


; e bs =0 


e a run 2 r :: C 


; e 8S >0 


a b run 3 e bs =0 


; e 85 <0 



CO 

Q. 



1 - 








0.4 0.6 0.8 1 

FIG. 2. Normalized density profiles p a (r)/p s of centers of 
the small spheres as a function of the distance r from the 
center of a single big sphere. Results are from computer sim- 
ulations for runs 1 — 3; i.e. ebs = 0, and e S3 is varied. The 
values of the e;j can be found in table |l| In this figure, as well 
as Figs. 3-6, the circles denote a repulsive e sa , the squares an 
attractive e ss and no symbol means no added e as . Solid lines 
denote no added e^, dotted lines et, B > 0, and dot-dashed 
lines e ba < 0. 
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FIG. 4. Normalized density profiles p a (r)/p a from com- 
puter simulations for runs 7 — 9; tbn < 0, e sa is varied. We 
also show the result for run 1, the pure HS case. 
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FIG. 5. 



Effective forces between the big particles, 



PFm, (r), shown for all 9 runs. The symbols are the same 
as in Figs. 2-4. The results are from direct computer simula- 
tions. Notice that the big-small interaction has the dominant 
effect. A repulsive results in a more repulsive ' (r) 
(dot-dashed lines) , while an attractive ebs results in a more at- 
tractive (r) (dotted lines). The small-small interactions 
have a relatively smaller effect. 




1 



r/c, 



1.3 



1.5 



FIG. 7. Comparison of theory to simulation for run 1: 
e bs = 0, e ss — 0. Simulation (solid line with symbols), su- 
perposition approximation (solid line), and DFT (dot-dashed 
line) results for the effective force and potential (see inset) as 
a function of r, the distance between the centers of the big 
particles, are shown. 
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FIG. 11. Same as Fig. (| but for run5: e bs > 0, e sa > 0. 
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FIG. 9. Same as Fig. 0, but for run3: e bs = 0, e ss < 0. 
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FIG. 12. Same as Fig. |, but for run6: e bs > 0, e sa < 0. 




FIG. 13. Same as Fig. but for run7: e bs < 0, e S3 = 0. 
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